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Abstract 

We consider a two-component regular cosmology bouncing from contraction to ex- 
pansion, where, in order to include both scalar fields and perfect fluids as particular 
cases, the dominant component is allowed to have an intrinsic isocurvature mode. We 
show that the spectrum of the growing mode of the Bardeen potential in the pre-bounce 
is never transferred to the dominant mode of the post-bounce. The latter acquires at 
most a dominant isocurvature component, depending on the relative properties of the 
two fluids. Our results imply that several claims in the literature need substantial 
revision. 



PACS: 98.80.-k, 98.65.Dx, 98.80.Es, 11.25.Wx 



Bouncing cosmologies have been proposed as possible alternatives to standard inflation 
in string-inspired (e.g. Pre-Big Bang pP, ekpyrotic/cyclic [2]) scenarios. In order to compare 
the predictions of these scenarios with observations, it is crucial to follow the evolution of 
cosmological perturbations from the initial vacuum-normalized state, through the bounce, 
all the way until decoupling. Nevertheless, there is no general agreement on the true influ- 
ence of a bounce on cosmological perturbations. In particular, it is not clear whether the 
Pre-Bounce growing mode of the Bardeen potential leaves any trace in the Post-Bounce, or 
whether it just completely matches a decaying mode Only in the former case can the 
claim that the ekpyrotic/cyclic scenario is in agreement with CMB data be defended. 

Studies of specific regular models apparently led to different results, suggesting that 
both possibilities can arise, depending on the specific model one assumes for the regular 
bounce. In spatially closed universes, a regular bounce from contraction to expansion can 
be triggered by the spatial curvature term in the Friedman equations [4 \. In this case, the 
transfer matrix relating Post-Bounce modes to Pre-Bounce perturbations depends on the 
momentum scale k, suggesting that non-trivial matching conditions hold. These models 
have the advantage of being completely embedded into General Relativity, but they need 
fine-tuning in the initial conditions in order to have a sufficiently long collapse. 

Alternatively, one can start from a spatially flat universe, avoiding the complications 
of spatial curvature. However, Friedman equations imply that a bounce from contraction 
to expansion necessarily violates the null energy condition (NEC). This violation can be 
induced by high energy and high curvature corrections to General Relativity or can be 
directly introduced by a ghost field dominating during the bounce. In this class of models, 
contradictory results have emerged up to now. Peter & Pinto-Neto jH] and Finelli [0] 
find some mixing between the growing and decay modes of the Bardeen potential, while 
Gasperini, Giovannini & Veneziano [7], Cartier [S], and Allen & Wands 9 find that the 
growing mode of the Bardeen potential in the Pre-Bounce exactly matches the decaying 
mode in the Post-Bounce. 

In this letter we report the results of the study of a wide class of regular bouncing 
cosmological models, containing, as particular cases, all the above-mentioned two-source 
regular models. We confirm that the Pre-Bounce Bardeen growing mode matches the Post- 
Bounce decaying mode. The only possible mode mixing is between the two components 
that are present in the cosmological description, as a consequence of the normal interplay 
between adiabatic and isocurvature perturbations. We present hereafter the key arguments 
and give the main results, postponing most details to a later publication |10j . 

We consider a two-component cosmology, where the first component dominates at early 
and late times, while the secondary one has negative energy density and triggers the bounce 
from contraction to expansion in a fully regular fashion. It is well known that ghosts or 
negative energy density fluids induce instabilities at a quantum level, but this has nothing 
to do with the evolution of classical perturbations, which is completely under control, as 
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we shall explicitly show. One may regard the secondary field as an effective term in our 
4-dimensional Friedman equations simulating corrections from high energy, high curvature, 
extra-dimensions, finite string size, or whatever else one can imagine. One can then look 
at the regular bounce obtained in such an effective description at a classical level using 
the paradigm of General Relativity, without worrying about quantum instabilities which 
should be properly addressed only within a yet unavailable complete theory of the bounce. 
These considerations confer a huge interest to ghost-induced-bounces as the only effective 
models of spatially-flat bouncing cosmologies which can be investigated within the frame 
and with the instruments of General Relativity. In some sense, they are the only models 
worth speaking about at the effective level! 

The background equations are 

m 2 = a 2 p (1) 
H 2 + 2H' = -a 2 p (2) 

p' + m{ P + p) = o, (3) 

where a prime denotes derivative w.r.t. conformal time rj, a is the scale factor, 7i = a' /a, 
and we have set 8irG = 1. The total energy density is p = p a + pb and the pressure is 
p = p a + Pb- The two fluids satisfy the continuity equation (jSJ) independently, i.e. we 
assume no interaction between the two fluids other than gravitational. We shall set rj = 
at the bounce, i.e. when TC changes sign. A typical background solution starts from a 
contraction driven by the first fluid, which ends when the secondary fluid's density becomes 
of the same order as that of the first fluid. We will call this time — % and normalize rj 
in such a way that r/b ~ 1. As we are restricting ourselves to spatially-flat cosmologies, 
the null energy condition (NEC, hereafter) is necessarily violated, i.e. p + p must change 
sign at some time f?NECi with — 1 < t^nec < 0. After the bounce the NEC is restored, the 
secondary fluid becomes subdominant again, and an expansion phase dominated by the first 
fluid follows. 

We shall consider only scalar perturbations, since these are both the most relevant and 
the most controversial. The perturbed line element takes the well-known form 

ds 2 = a 2 (r]) {(1 + 2(j))dr] 2 - IB^dx 1 

- [(1 - 2V>)% + 2E jij ] dx i dx j } , (4) 

while the perturbed energy-momentum tensor of the first source reads: 

frpy = ( Pa + Spa ~{p a + Pa)V a ,i \ ^ 

V \(Pa+Pa)Va,i -{Pa + $Pa)$ij /' 

where V a is its velocity potential. The same expression holds for the secondary fluid, the 
subscript a being replaced by b. 
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In order to keep our treatment as general as possible, we introduce the following gauge- 
invariant relations among the perturbations of the sources: 

bp a = c 2 Jp a + a(p a + p a )V a ; a = a a a ; Sp b = (? h bp h . (6) 

Pa + Pa 

Notice that if c 2 ^ p' a / p' a , the additional dependence on the velocity potential generates an 
isocurvature mode 

T a 5S a = 5p a - {p' a /p' a )5p a = a{p a + p a ) (V a ~ 5pa/p' a ) • ( 7 ) 

It is a trivial exercise to prove that a scalar field <p with a self-interaction potential V(ip) 
satisfies these relations with c 2 = 1 and a = — 2a 2 V t p/ip' . Therefore, these relations include 
perfect fluids and scalar fields as particular cases of a general class of sources and ensure 
a wide generality to our treatment. Allen & Wands' model [5] is recovered with the scalar 
field choice for the first fluid and c 2 = 1. The perfect-fluid model studied by Finelli |H| is 
recovered setting a = 0, = 2(? a + 1/3. The model by Peter & Pinto-Neto [H] corresponds 
to a = 0, c 2 = 1/3, = 1. The generalization to the case where an intrinsic isocurvature 
mode for the secondary fluid is present poses no major difficulties, but will not be considered 
here for simplicity. 

The first subtlety in the study of perturbations in cosmological bounces is the choice 
of a gauge where all variables stay finite and small enough for the linear theory to apply. 
In particular, we have to be sure that the perturbative variables we are considering will 
not diverge at any point. From the physical point of view the requirement is that all 
components of the metric and of the energy-momentum tensor stay finite at least in one 
convenient gauge. If we assume that such a gauge exists (if it does not then the model 
under consideration cannot be a stable solution of Einstein equations), we deduce that 4>, 
tp, B, E, dpi, 5pi and (pi + pijVi stay finite and small in that gauge. We will refer to it 
as the regular gauge. Note that the total velocity potential V is allowed to diverge in the 
regular gauge at rj = t/nec when p + p = (NEC violation), since what matters is that the 
energy-momentum tensor components stay finite. Our attitude will be just to assume that 
such a regular gauge exists and to build gauge-invariant quantities out of perturbations in 
this gauge. 

The Bardeen potential \P = ip + 7i(E' — B) is indeed a first regular combination. The 
widely used £ = ip + TtV may diverge at t/nec- However, we can replace it by 

C = (H 2 -H')(i> + HV), (8) 

which is guaranteed to be regular, since 2(H 2 - H')V = a 2 (p + p)V is finite in the regular 
gauge. We then introduce individual gauge-invariant variables 

Ca = \a 2 {Pa+Pa){i>+HV a ) (9) 

H> a = ^a 2 (p' a if> + H5p a ), (10) 
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and similarly for the b fluid. All these variables are regular by construction. From the 
definition of £ and the 00 component of the Einstein equations, we have 

C = Ca + 0>; (3C + kH)n + * a + v b = o. (n) 

The perturbations of the continuity equations of each fluid can be combined to give 
coupled second order equations for ( a and £&■ Introducing the canonical (Sasaki-Mukhanov) 
variables \H\ v a = ( a /[c a (Pa + Pa) 1/2 W] and v b = ( b /[c b (p b + p b ) 1/2/ H], we find 

v" a + (c 2 a k 2 - 4/z a ) v a = 0{p b / Pa ) 1 / 2 (12) 

V'l + (elk 2 - 4/ Zb ) V b = Oiph/pa) 1 / 2 , (13) 

with z a = a?{p a + Pa) 1 / 2 1 (c a 7~C) and z b = a^ 1+3c b\p b + p b ) l l 2 jc b . These equations also 
identify c 2 a and as the correct definitions of the sound speed for the two components. 

These equations decouple only in the limit p b /p a — ¥ 0. Both Mukhanov variables are 
normalized to vacuum fluctuations and thus are of order l/\^2k in the asymptotic past. 



However, we have ( b = \/(p b +p b )/{pa + PajQa and a hierarchy can therefore be established 
between these two variables in this limit. In this sense, the 0(p b / p a ) l l 2 contains all terms 
that are irrelevant in the asymptotic past, including couplings between the two variables. 

Although it is not necessary, let us assume for simplicity that the Pre-Bounce can be 
described by a power law a(r/) ~ (—rj) q . Defining T = ^j, the following general relations are 
implied by the background equations 

- 2 a ~*( 1+3c S- 2 . (14) 



* 1 + 3T' -rj 

If a = 0, we get F = c 2 ., which is the case of a perfect fluid. A well-known case of a source 
where a contributes to the background is the case of a scalar field with an exponential 
potential e~ Xip . Then c 2 . = 1, but T can assume any value depending on A. This case is 
discussed in Ref. f9.. A stiffer 77-dependence for a would make it rapidly subdominant, 
while a softer one would render the vacuum normalization problematic and will not be 
considered. 

The secondary field evolves according to its continuity equation p b ~ (— rj)~ 3q ( 1+c b\ 
Exploiting the fact that the two fluids have comparable densities at time rj = — 1, at earlier 
times we find p b /p a — (—r])~ q ( 1+3c t>") +2 . 

In the power law regime, Eqs. (|12j) and give 

^\C a H${caVk\n\) (is) 
v b = ^\C b H${c h Vk\ n \). (16) 

with u a = \ — q, u b = ^(3c^<7 — 1 — q) and is the Hankel function of the first kind; C a 
and C b are pure numbers of order 1. The initial conditions for £ a , ( b , fy a and ^ b can be 
derived from those we just imposed on Mukhanov's variables. 
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From the perturbations of the continuity equations, it is straightforward to write a set of 
first-order differential equations for f, ^, fy b . However, it is more useful for our purposes 
to write these equations in their integral form 



a 2 



a 2 - 



ci(fe) , / , 



(17) 



C = a 2 (p a +Pa) 



C2(fc) 



a 2 (/Oa +Pa 



-k 2 ^di] 



a 



a 2 (p a + Pa 



<bdrj + 



a 2 (p a +Pa 



(18) 



n 



-Cbdi] 



+ 



- 3+3cl (p b +pbl k ^ d71 



27i 

Cb =a 1+ ^H( Pb + p h ) 



(19) 



c 4 (k) + 



a 1+3c tH(p b +p b 



■^bdv 



+ 



1 — Sr 2 
a 1 dC 6 ~ 



2H 2 



Qdr] 



(20) 



Now we are ready to analyse the evolution of the modes that are well outside the horizon 
at the bounce (k <C 1). The discussion can be made in terms of the two parameters T and 
Cj, which express the background evolution of the two fluids. Not all values of these two 
parameters are relevant to bouncing cosmologies. Indeed, as a first constraint, the condition 
that Pb/pa must grow during contraction (otherwise there is no bounce), excludes the half- 
plane c 2 < r. A correct vacuum normalization requires c 2 ,c 2 > 0. Finally, we require 
r > —1/3, in order to discard super inflationary backgrounds. For the sake of generality we 
do not impose the causality constraint c 2 , c 2 < 1. In Fig. ^ we show the allowed regions in 
the cases a^O and a = 0, the forbidden ones being shaded. 

The four constants C{{k) can be determined by matching them to the long- wavelength 
limit of Eqs. (1T51) . (fTK|) : 



ci{k) ~ k Va , Qj(Jfe) ~ k~ Ua , c 3 (k) ~ k~ v \ Ci (k) ~ k 



(21) 



Now we have to compute the evolution of the individual terms in Eqs. (|17JI -I|2U |) and 
find the dominant contribution at horizon re-entry after the bounce. We will discard any 
numerical factors of order 1, just keeping the dependences on k and 77, the only dimensional 
parameters in our analysis. In these equations, each variable is expressed as the sum of a 
"homogeneous" mode and several integrals, one for each coupling term. The homogeneous 
mode of each variable persists through the bounce and after the bounce, and thus poses no 
difficulty. To evaluate the integrals, we note that, far from the bounce, we can replace all 
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Figure 1: The parameter space in the general case (left) and in the perfect fluid case (right). 

the functions by their asymptotic power law behaviours . Thus each integrand /(??) assumes 
the generic form Ci(k)\rj\ s and, for each integral, we generically have, in the Pre-Bounce: 



Ci(k)\rj\ s dr] ~ Ci(k)\rj 



s+1 



rj < -1, 



(22) 



where the integration constant is set to zero in order to match, in the far past, the asymptotic 
solutions implied by (|T5|) . ffTTTf) . 

To evaluate the integrals in the Post-Bounce, we can split the integration domain in 
three intervals, covering the full Pre-Bounce, the Bounce and the Post-Bounce separately. 
We have 

-l l 
Ci{k)\ri\ s dr} ~ Ci(fc), / f(r])dr}~Ci{k), 



-i 



d(k)\r]\ s dij ~ a(k) I — 77 



(23) 



The first integral is just Eq. (|22|) evaluated at 77 = —1, while the third is evaluated using the 
power-law expansion for the integrand 7(77) in the Post-Bounce. More subtle is the bounce 
contribution. Since all the curvature and time scales are fixed by the time normalization 
choice rjb = 1, the only scales that can be introduced by this contribution are k and 77. 
Yet, the integrand is a combination of k- independent backgrounds, and integration over 
the finite bounce interval eliminates 77 from the final result. Hence the bounce contribution 
is of the same order as Ci(k). An apparent complication comes from the divergence of 
some integrands in Eqs. (|17|) - ([20j) at the bounce. Of course, these divergences do not 
affect the regularity of the gauge-invariant variables, since all the integrals are multiplied 
by proper factors, which compensate these divergences. Mathematically we can further split 



6 



the bounce integrals in two pieces, before and after rj = 0. The limiting value of the gauge- 
invariant variable at 0~ fixes the initial condition for the second integral starting from + 
(recall that our integral equations are just another representation of differential equations, 
which admit a regular expansion for the solution at r/ = 0). Afterwards, the dimensional 
argument can be applied safely. 

At late times (77 S> 1) the dominant term depends on the value of s 



Calculating all the integrals recursively, we can build a ^-expansion for each of the four 
independent modes Ci(k). Evaluating the solutions in this way obtained at the horizon 
re-entry 77 ~ 1/k, we can pick out the dominant contributions and read their spectra. 

Let us first discuss the general case where a is different from zero. This includes bounces 
with tracking scalar fields and all asymmetric bounces. We can distinguish two regions in 
the plane (I\cf), illustrated in the left panel of Fig. ^ In region A (defined by c 2 < 1), ^ 
and C develop a typical adiabatic spectrum, determined by the dominant fluid. We have 
\I> ~ C ~ c\{k). In region B, covering the rest of the plane, and Q acquire the spectrum 
of the secondary component, so that ty B ~ Q B ~ c%{k). The individual variables Q a and 
Cb follow the behaviour of £. Notice that on the border line (c 2 = 1), the two spectra 
coincide. It was this line that was explored in Ref. In the particular case a = 0, we 
have r = c 2 and the bounce is always perfectly symmetric. This includes bounces with two 
perfect fluids or scalar fields dominated by their kinetic energy (compare Refs. [HIE])- The 
time-reversal symmetry induces several cancellations in the integrals, allowing the creation 
of an intermediate region C, defined by (2 + c 2 a )/3 < c\ < mm(5/3 + 2c 2 ., 2 - c 2 a ). In this 
region, £ a and Cb keep their original spectrum, without mixing. Nevertheless, $ and £ are 
dominated by the first component, so that \P ~ £ ~ c\{k). Finally, in region B, ^ and £ 
get an extra k 2 factor w.r.t. the case 

In no case does the spectrum k~ 2 c\{k) ~ k~^~ q of the Pre-Bounce growing mode of 
\I/ survive at the Post-Bounce horizon re-entry. All the above results have been confirmed 
by numerical integration of the equations ^U]. Our results support the conclusion that a 
smooth bounce cannot generate a scale-invariant spectrum via the mode-mixing mechanism 
advocated in One can argue, of course, that the way string theory will actually describe 
the bounce has no resemblance to a regular bounce described by General Relativity. In 
that case, understanding the spectrum of perturbations in such models will have to wait 
until a full stringy treatment of the problem is found. On the other hand, claims that 
such a treatment will automatically lead to a flat spectrum through mode-mixing appear 
unjustified. 

Some authors, working within specific models, have reached conclusions that appear 
to contradict our general analysis. In Ref. 5 j, Eq. (26) provides an analytical solution 




s < -1 
s > -1. 



(24) 
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fiv) = f?/(l + f/ 2 ) 2 for perturbations outside the horizon which is valid both far from and 
through the bounce. However, the authors use an approximate form (f ap p(i]) = I/ 7 ? 3 ) of 
this solution far from the bounce and then match f app to / just before the bounce. Had 
they used / from the beginning throughout the bounce, they would have found no mixing. 
No surprise, therefore, that the mixing they find is at order rj~ 5 , which is just the next to 
leading order term in f(rj), which is missing in f a pp(v)- Had they kept the r/ -5 term in 
fappiv) they would have found mixing at the order rj~ 7 . We thus conclude that the mixing 
they find is just an artifact of their matching procedure. 

In the case of Ref. jH], the limit 7 —* {k —* in our notation), used to evaluate the 
final spectrum, just selects the infrared modes at the bounce. However, the final spectrum 
should be evaluated at the horizon re-entry, after the decaying mode has gone away; it 
would correspond to 7 — > while keeping 72; (our krj) of order 1. Finally, a recent pa- 
per |12| studying the slow contraction of the ekpyrotic model in the synchronous gauge 
agrees with our general conclusions. However, these authors ignore a possible divergence 
of the £ variable at the NEC-violation point, and this appears to restrict significantly the 
applicability of their argument. 

V.B. thanks the CERN Theory Unit for hospitality. We acknowledge useful discussions 
with F. Finelli, J. Hwang, V. Mukhanov, P. Peter and D. Schwarz. 
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